Exact solutions of a class of = 1 quantum Ising spin models 
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We propose a hole decomposition scheme to exactly solve a class of spin-1 quantum Ising models 
with transverse or longitudinal single-ion anisotropy. In this scheme, the spin-1 model is mapped 
onto a family of the S = 1/2 transverse Ising models, characterized by the total number of holes. A 
recursion formula is derived for the partition function based on the reduced 5=1/2 Ising model. 
This simplifies greatly the summation over all the hole configurations. It allows the thermodynamic 
quantities to be rigorously determined in the thermodynamic limit. The ground state phase diagram 
is determined for both the uniform and dimerized spin chains. The corresponding thermodynamic 
properties are calculated and discussed. 
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I. INTRODUCTION 

The phase transition driven by quantum fluctuations 
is one of the fundamental issues in quantum many body 
systems. A number of novel phenomena associated with 
the transition such as the quantum critical behavior have 
been observed in a variety of condensed matter materials 
P, 0] ■ One of the prototype model systems exhibiting the 
quantum phase transition is the one-dimensional spin- 1/2 
Ising lattice with a transverse field, namely the transverse 
Ising model (TIM) a, 4, 1] , defined by 



H' 



TIM 
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(1) 



jj is the spin operator at the site j on a one- 



where Sj 

dimensional lattice of length L. The transverse field 
h introduces quantum fluctuation to the system, lead- 
ing to a quantum phase transition from the ferromag- 
netic/antiferromagnetic ordered to the paramagnetic dis- 
ordered states above a critical value he = J/2. Actually 
the model is equivalent to a free spinless fermion system 
and can be exactly solved by applying the Jordan- Wigncr 
transformation!^ 0|. Based on the exact solution, all 
physical quantities, including the ground state energy, 
the low-energy excitations, the specific heat and other 
thermodynamic functions can be evaluated. This pro- 
vides a thorough understanding of the quantum critical 
behavior of this system. 

However, in real materials, the moments of atoms may 
be larger than 1/2. Despite of immense efforts in the 
past three decades it is still very difficult to find exact 
solutions of the S" = 1 or other higher spin quantum 
Ising systems. This is partly due to the existence of the 
spin neutral states S" = (called as holes hereafter) in 
addition to the two spin polarized states = ±1 at 
each site in the 5* = 1 spin chain. In these models, a 
hole can decay into a pair of polarized spin states, and 
vice verse, thus rendering the Jordan- Winger approach 
invalid in exactly solving the S — 1 TIM. 

In this paper, we study a class of one-dimensional 



S = 1 quantum Ising model, defined by the following 
Hamiltonian 



H=-J2iJ^S3Shl+f3) 



(2) 



Where, fj ^ D^{S^y + D'^{S}I)^ + D^{S^)^ is the single- 
ion anisotropy term with site-dependent D'j{a — x, y, z). 

The above model has a classical limit where Dj = Dj. 
This corresponds to the Blume-Capel modelQ. So the 
present model can be regarded as the quantum general- 
ization of the Blume-Capel model. The simplest quan- 
tum case is the uniform chain defined by[i| 



QIM 



E 



(3) 



In two or higher dimensions, this kind of quantum 
Ising model with single-ion anisotropy was studied by 
a number of authors, based mainly on the mean-field 
approximations [lol [llj. In particular, the ground state 
of the model ^ was shown to be equivalent to the 
S" 1/2 TIM defined by Eq. (P) p^. Such equivalence 
is also valid for quantum Ising models with bond- and 
site-alternationsil3| or geometrical frustrations, such as 
a fully frustrated spin-1 Ising Delta-chain[14|]. 

The purpose of the present paper is to study the physi- 
cal properties of the model described by Eq. ([2]) based on 
the exact solution. The key idea in solving the proposed 
S — 1 model is to divide the total Hilbert space of the 
S = 1 system into a number of subspaces labeled by the 
number of holes. This is what we call the hole decom- 
position scheme (HDS). This HDS was developed in our 
recent work [9| where a recursion approach based on the 
HDS is suggested for the uniform chain. In the present 
paper, we shall give more comprehensive investigations 
for various properties of the 5* = 1 model, including the 
case with dimerization [Tsj . In particular, we show that 
for a given hole configuration, each sub-lattice system 
with purely polarized spins can be exactly solved in that 
case. Based on the exact solution we study how the quan- 
tum phase transitions and thermodynamic properties are 
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affected by tfie interplay between the dimerization and 
the single-ion anisotropy. Depending on the strength of 
dimerization, we find that the system undergoes a quan- 
tum phase transition where the criticality is the same as 
that of the uniform 5 = 1/2 TIM. 

We note that the spin- 1/2 TIM can be realized in cer- 
tain low-dimensional magnetic materials p^ . [iTj . For the 
systems with local moments larger than 1/2, the single- 
ion anisotropy generated by crystal fields and the dimer- 
ization may also become important [l^. The 5=1 TIM 
with the crystal field splitting was used to describe the 
ferroeletric transition in SrTiOsflo'l. In a class of quasi- 
one dimensional spin chains, such as [Ni2(Medpt)2(i^- 
ox)(H20)2] (0104)2 •H2O, the magnetic Ni2+ ion shows 
not only a single-ion anisotropy [20, [2l|, but also a, S — 
1 bond- alternating pattern, where Medpt is the bis(3- 
aminopropyl) methvlamine [2^ [23j . Recently, cold atoms 
or polar molecules in optic lattices were shown to be 
ideal systems to realize various quantum spin models [2^. 
In particular, the spin-1 models can be implemented 
by trapping polar molecules where the spin degrees of 
freedom can be described by the hyperfine vibrational 
states j25|. 

This paper is organized as follows. In Sec|TT] and 
Sec lIIIi we discuss some general properties of the model 
and introduce the HDS. In Sec lIVl we solve exactly the 
Hamiltonian in the presence of dimerization. In SccfVl 
we study the low energy excitation spectra and quantum 
phase transitions based on the exact solutions. In Sec lVIi 
we discuss in detail the recursion method introduced in 
Ref. for evaluating thermodynamic quantities. Fi- 
nally, we give a summary in Sec. IVIII 

II. HOLE DECOMPOSITION SCHEME 

Let us consider a S=l Ising lattice with single- ion 
anisotropy, defined by Eq. {Sj,Sj,Sj) are the 

5 = 1 spin operators at lattice site j = 1,2,...,L, 
with the lattice length L. The uniform classical Blume- 
Capel model [8] corresponds to the symmetric case with 
— _D| = D. Because only two of these terms 
are independent, we shall mainly consider the quantum 
case with Dj = Dj and = D| = 0, without losing 
generality. The case with D| ^ will be discussed later. 
We shall mainly focus on the dimerization case where 

J2j-1 = Jl, J2j = J2, D2J-I = £>!, D2j = D2. (4) 

By definition, one has {S]f = \{S+ S+ + SJ SJ + 
Sj~S~ + S~Sj~). So it is straightforward to show that 
(5J)^ does not couple 5| = state to 5| = ±1 states. 
Thus (5J)2 commutes with (5|)2. This leads to the fol- 
lowing theorem: 

Theorem 1. When 5 = 1, the total hole number opera- 
tor Nq — i — commutes with the Hamiltonian 
for any site- dependent Jj and Dj 

[No,H]=0. (5) 



It means that the total number of holes is a conserved 
quantity if 5 = 1. Consequently, the eigenstates of H can 
be classified by the eigenvalue p oi Nq. In the discussion 
below, we will call the subsystem with p holes as the p-th 
sector. 

Let Hp be the Hilbert space of the p-th sector, the 
total Hilbert space is then given by the sum of all the 
subspaces: 

n = nQ®ni®---®HL. (6) 

A complete set of the eigenstates in the p-th. sector forms 
a sub-band of the whole spectrum. Generally the lowest 
eigen-energy of the p-sector satisfies the following theo- 
rem: 

Theorem 2. Let E{p, L) be the eigen-energy, corre- 
sponding to the eigenstate |p, i) of H, and Eq{p, L) be 
the lowest eigen-energy in the p-th sector, then the fol- 
lowing inequality holds 

E„ip,L)<Eoip+l,L). (7) 

This relationship is an extension of the Lieb-Mattis 
Theorem derived initially for the uniform systems with 
ferromagnetic Ising couplings [26*1. It holds still no mat- 
ter whether the Ising couplings arc ferromagnetic or an- 
tiferromagnetic. It indicates that in the absence of the 
_D|-term the ground state always lies in the p = sector 
and the energy spectrum of the system has a hierarchy 
structure. 

There are two kinds of excitations in the system. One 
is the fermionic excitation within a given sector. The 
corresponding excitation energy is defined by E(p, L) — 
Eq {p, L) . The other is the hole excitation and the exci- 
tation energy with respect to the ground state is given 
by E{p,L) — E(){0,L). There are two kinds of minimal 
excitation gaps corresponding to these excitations 

A*,^) = E,ip,L)-Eoip,L), (8) 
Ajf) = Eoip,L)-EoiO,L). 

Fig. [1] shows schematically the hierarchy band structure 
of the system. Within each bands (or each boxes shown in 
Fig. [T]) there are fermionic excitations, with the minimal 

gaps Aq^-* . While the minimal gaps of the hole excitations 
aJj''^ increase with p when the longitudinal anisotropy 
= 0. 

HI. MAPPING ONTO THE SPIN-1/2 TIM'S 

The eigenstates of the Hamiltonian Eq. ^ 
can be generally expressed as |^) = 

X^mj -^™l.™2,--- "^2, • • • TOl), with Fmi,m2,--- ,mL 

being the wave function. As the total spin and its 
z-component are not conserved, the summation runs 
over all nij = 0, ±1 states. However, by Theorem 1, 
the local holes are good quantum numbers and can be 
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FIG. 1: (Color Online) The schematic picture of the energy 
band hierarchy. The boxes represent the sub-bands of the 
corresponding p-ih sectors. Aq''^ and A^''^ are the minimal 
gaps for the fermoinic and hole excitations respectively in the 
p-sector. Several explicit results of these gaps will be given in 
Section IV. 



regarded as non-magnetic local impurities embedded 
in the S — 1/2 Ising system. Assuming in the p-th 
sector the holes are located at x„ {n — 1,2, ■ ■ ■ ,p), 
then the corresponding eigenstates can be expressed 



Yl.^i<^\mj), where 



±1 if j ^ Xn, and Tfi^i 



0. 



In the p = sector, there are only two spin states 
at each site, corresponding to mj = ±1 respectively. 
They are in one-to-one correspondence with the states of 
the spin- 1/2 (Pauh) operators tT|: 5'||mj) <^ (T||mj) = 
rhArhj). One has then the mapping relationship: 
Sfsf ^ 2af, SfSj - 1 ^ ±cr|. Thus, (S'J)^ acts 
like (1 -I- o-j)/"^- The original Hamiltonian, when acting 
on the p — subspace, has the following reduced form 



(9) 



This is just the spin-1/2 Ising model with bond (site)- 
dependent Ising couplings and transverse fields. 

Now let us turn to the p = 1 sector. If 
the hole is located at the site xi, the cor- 
responding state can be written as |^'(i)) = 

Ea, F^■%^,-^.,o,rh^,+^,■■■\r^^l ■ ■ ■ m^^^-iO^.m^^+i ■■■mL). 
Now because the bonds connecting the hole are broken, 
the reduced Hamiltonian, obtained by acting the original 
one on the p = 1 sector, is given by 



(10) 



Where i/^'^ is the Hamiltonian of the spin-1/2 TIM seg- 
ment of lengths I imposed by the open boundary condi- 
tion. 111-0 (xi — l)) and \tp (L—xi)) are the eigenstates of 
the segments ff^'g and i^_r^^y respectively, then 

= - 1)) ® |0,,) ® |0(L - X,)). (11) 



is the eigenstate of . This HDS can be easily gener- 

alized to the multi-hole sectors. For instance, the reduced 
Hamiltonian in a p-hole sector is given by 



iP,L) 



H{Q,xi-l} +-^(0,1:2-2:1-1) 



(12) 



where the positions of holes. Similarly, the eigen- 

state of H(^p ^ can then be expressed in terms of those 
of the individual segments. 

Therefore, all the eigenstates of the 5 = 1 QIM ^ can 
be obtained by solving a set of the spin-1/2 TIMs. For 
the p-th sector, the number of the S=l/2 TIM segments 
are p + 1 or p, depending on whether the original chain 
is periodic or open. Note that the hole positions may 
vary along the chain, so there are many different hole 
configurations in a given sector. The decomposition of 
the total Hilbert space into the sum of subspaces can be 
formally represented by [2 © l]*"-^ = 2'^^ 2®(^-i' (g) 1 ® 
2®(L-2)^ -1^(812^. . .Qi'S'i^ where the dimension of the S = 

1 system is given by dimH = (2 + 1)^ = Ylp^o 2^"^C^ = 

IV. DIMERIZED SPIN CHAIN 

A. The non-hole sector 

In this section, we follow the standard approach intro- 
duced in Refs. [a, [23| to diagonahze the p = sector of 
the S=l model, Eq.(|9]). We introduce the fermion opera- 
tors Cj and Cj and use the Jordan- Wigner transformation 
to rewrite Eq.Q as follows 



(OX) 



E 



CjBjqCq) 



+HpB, 



(13) 

where Hpb = Jl{c^lc{ + c\ci - clc\ - clCi){Kl+i + 1), 
JjSj+i,q ^ Jj-i^j~i,qi a.nd Bjq = 



A ■ — 4 ■ — 



-B„ 



'Dj5j,q 



3"3+i,q ■ 



Jj-iS 



Notice that Hps is the 



boundary term (which disappears for the open chains) 
and can be neglected in the thermodynamics limit. 

Next, we introduce the Bogoliubov transformation as 
in the following: 



Vk 



kj Cj 



hkjC 



Y.i9lA + ^k,c,). (14) 



Where, rjk and 77]^ are the fermionic quasi-particle oper- 
ators with quasi-momentum k. They satisfy the usual 
anti-commutation relations. While, gkj and hkj are co- 
efficient matrices, which should be complex in general. 



Then, i?(o,L) becomes 



H{o,L) = E ^(f'^ivk + 1 E 

fc 



k 



(15) 
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If we define ($fe)j = gkj + hkj, (*fc)j = gkj - hkj, the 
eigenvalue A(fc) can be solved by 

M,,{^k)g = A'{k){<t>k)r,M,,{^k)g = A2(fc)(vI/,.)„ (16) 

with M being a symmetric matrix defined by M = 
{A - B){A + B), or M^, = {DjDq + 4J,_i J,_i)(5,- , + 
2DqJj-iSj^i^q + 2DjJg-iSj^q-i. For the dimerized (or 
alternating) chain, the M-matrix depends on the bound- 
ary conditions and the parity of the chain length (even or 
odd) . For simplicity, we here consider the case of periodic 
chain with even L, where 



M 



( ai 


bi 


• 


• 


b2 




02 


62 • 


• 








b2 


ai ■ 


• 











• 


• ai 


61 


\h2 





• 


■ bi 


a2 



(17) 



with ai ^ Df + 4J|, bi = 2DiJi, a2 ^ Dj + AJl, and 
62 = 21)2 J2. 

In order to solve Eg. pop with this M-matrix, we em- 
ploy the following Ansatz 



(*fe)2j 
(*fe)2j + l 



Ae exp[ifc(2j)], 
Ao exp[i/c(2j + 1)]. 



(18) 



The ratio r = Aq/ A^ is then a measure of dimerization, 
determined by 



ai — a2±W 



(19) 



where 



W EE v/(a2 - ai)2 + 4(&i - &2)2 + I66162 cos2 fc. (20) 

Note that when ai = 02 and 61 =62, we have r = ±1, as 

it should be for the uniform chain. 

The eigenvalue A'^{k) is then obtained by A'^{k) = 
± so that one has four branches of quasi-particle 

excitations: 



A±rik) 



±r 



with r 



1,2, r2 



1 - (-l)Vl -^ + '9cos2A: , (21) 



pp4 



DjD^ + 16Jf J|, gr^ = 8D1D2J1J2, and k = 2nm/L, 
{m = -L/2,---L/2- 1). 

It should be noticed that A±r(k) (for r — 1,2) are 
invariant under the shift k k + tt. Thus the quasi- 
momentum k is constrained in one of the following 
regimes: (1) k e [-7r/2,7r/2) or (2) k € [-n,-n/2)U 
[7r/2,7r). Here, we choose k to be in the first regime 
for Aj-i(fc) and in the second regime for A±2{k). Fig. [2] 
shows the typical energy dispersions for a dimerized sys- 
tem with Q > 0. The case with Q < can be obtained 
by a reflection under the shift. 



< 





-~^,(k) 







-ji/2 
k 



7l/2 



FIG. 2: The representative fermionic excitation spectra in the 
non-hole sector of the 5" = 1 QIM with Q>0, Ji = 1.01, J2 = 
1, D\ = 1.82, D2 = 0.9. This kind of excitations for a sector 
with fixed hole number is similar though the quasi-momentum 
is determined by the corresponding secular equation. 



B. The single- hole sector 

In the p — \ sector, the hole breaks two bonds con- 
necting to it. For convenience, we assume that the hole 
is located at the site L. The reduced Hamiltonian can be 
rewritten as 



^(1 



(O.L-l) 



(22) 



where i/^'^ ^ 



is defined by Eq. but with the open 
boundary conditions. The procedure for diagonalizing 
i/^'p is the same as for H(^Q ]^y The eigen wavcfunc- 
tion is determined by Eg. pop in a similar way. The only 
difference is that now the M-matrix is a (L — 1) x (L — 1) 
matrix, with Af^-i 1 = Mi i_i = 0, and Mi 1 = ao(= 
Dl). 

Due to the open boundary condition, the Ansatz 
Eq. (US]) is no-longer valid. Intuitively, a reflection wave 
(g-jfcj) Tffiw be induced at the boundary in addition to the 
incoming wave (e''^-' ). Thus we suggest another Ansatz 
for the dimerized open chain as follows: 



(*fe)2j 
(*fe)2jH 



^2ikj 



tee 



-2ikj 



^A_^ 



(23) 



With this Ansatz, we find that t and A{k) take the same 
form as those defined in the p = case. The reflection 
coefficients, to and tg, are given by 



te = 



iLk ' 



tn 



{hie"-^ + b2e^^^)e-'-^^' 

Moreover, the quasi-momentum k is now determined by 
the following secular equation 

262 [bl + bl + 2bib2Cos2k\ sin Lk 

= (ao — ai) [61 sin(_L — 2)fc + 62 sin Lk] 

X [(ai - 02) ± W] . (25) 
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The equation is symmetric under k —k, thus we only the ground state is given by 
need to solve the complex or positive /c's. The solution 
can be further simplified if 61 = 62 ( or DiJi = D2J2) 



where to ^te. Ao ^ \ - JT^ ~ St, (27) 



C. The multi-hole sectors 

The previous approach is extended to the subsystems 
or sectors with more holes. The holes break the Ising 
couplings, leading to a series of segments of the spin- 
1/2 TIM's. For the periodic chain, the number of these 
segments is equal to the number of holes p (including 
the segment of zero length). If the holes are located at 
(xi, X2, • • • ,Xp), then the corresponding reduced Hamil- 
tonian is given by Eq.®, i.e., if(p,L) = J2n=i H {0,1^)- 
Y7n=i ^^r.^ with ln = Xn- Xn-i - l(here, xq = Xp - L) 
being the length of the n-th segment and 

i=a:„_i + l 

j=2;„_i + l j=2;„_i + l 

Each segment Hamiltonian iJ'(0,^„) can be diagonal- 

ized as H'iO, Q = Efc A(fc) {^rjU^ - k) "EJ^i.+i ^, 
where k is the quasi-momentum satisfying the secular 
equation Eq. (pS)) (by replacing L by In)- Notice that 
associated with the fixed length of the segment there 
are four different kinds of configurations, depending on 
whether the two edge holes are located at odd or even 
sites. 



V. THE LOW ENERGY SPECTRA 
A. Fermionic excitations in the non-hole sector 

In the previous section, we show how to diagonalize 
the reduced Hamiltonians of different sectors. In a given 
hole sector, there are four branches of quasi-particle ex- 
citations, given by Eq. (PT|) . Obviously, the two negative 
branches A-r{k) (r — 1,2) will be filled in the ground 
state of that sector. 

In the p = sector, the ground state is given by 

IV'o) = Ukk''nlukvU,k'\^)^ where 77^^ (r = 1,2) are 
the fermionic quasi-particle operators in the branches 
A-r{k), k and k' are the allowed momenta for r = 1,2, 
respectively. The ground state energy is given by Eo{p = 
0,i) = ^EfcA_i(fc) + iEfe'A_2(fc'). 

The low energy excitations can be obtained by apply- 
ing operators r/^ j., or rj^2,k' on the ground state. The 
excitation energy is given by A2(fc')(= — A_2(fc'))- When 
Q > 0, the energy gap between the lowest excitation and 



where 6± = ^/\D^D^±Ij^\. If \Di\ = |i?2| = 1^1 
and |Ji| = IJ2I — \J\, one has Aq = 2|J|||^ — 1| with 
A = 2J / D, reproducing the result of the uniform chain. 
In this case, the branches 1 and 2 connect smoothly. In 
the presence of dimerization, however, the two branches 
will split and produce a dimerization gap 




B. Single hole excitation 

The ground state of the p — I sector has L-fold degen- 
eracy in the uniform chain, because the energy does not 
depend on the position of the hole. But in the presence of 
dimerization the degeneracy will be L/2-fold. While the 
momentum-dependence of the spectra can be determined 
as in the p = sector, the quasi-momentum k must sat- 
isfy the secular equation associated with the open chain 
of length L ~ 1. 

In the lowest state of the p — 1 sector, the bands 
with negative energies are fully filled. The correspond- 
ing energy is then given by Eq{1,L) |EfcA-i(fc) -|- 
^j., A_2(fc')] — i max (£>i, £)2). The last term is con- 
tributed from the hole which may locate at either even 
or odd sites. 

According to Theorem 2, the ground state of the orig- 
inal Hamiltonian should be the one in the p = sector. 
So in addition to the fermionic excitations, inclusion of 
the p ~ 1 sector (adding a hole to the system) will induce 
the hole excitations. The minimal hole excitation gap is 
given by A|^'^ = = Eq{1,L) - Eo{0,L). The typical 
behavior of A^ (as a function of D2) is shown in Fig. [3] 
for fixed Ji = 10, Di = 2.0, and L = 2000, with either 
periodic or open boundary conditions. 

In order to understand the difference in for the 
periodic and open chains, we consider a special case 
where the chain is uniform and at the critical point, i.e., 
Di ^ D2 ^ D, Ji = J2 = J, and A = 2J / D = 1. 
In this simple case, the energy spectrum is given by 
K{k) = -|D|VTTA2T2A cosfc, and the allowed quasi- 
particle momenta are k = 2™:^ ^ [rn — — L/2, • • • , L/2— 1) 
for the period boundary condition, and k = yfq^, 
(to = 1, • ■ • ,L) for the open boundary condition. Thus 
for sufficiently large L we obtain Ah ~ 0.1361? and 
Ah ~ 0.3181? for the open and periodic chains, respec- 
tively. Their difference, — 0.182L), does not change for 
larger L. So it is the surface energy cost to turn into an 
open chain. 
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FIG. 3: The hole excitation gap Aj, as a function of D2 for 
the system with periodic (left panel) or open (right panel) 
boundary conditions. Three curves corresponding to J2 — 
0.5, 1.0 and 1.5 are plotted and Ji — 10 and Di = 2.0 are 
fixed. 



FIG. 4: (Color Online) Ground state phase diagram of the 
S — 1 QIM. The f2-axis and Ti-a^xis describe the dimerization 
strength and the competition parameter, see Eq. (|30|I in the 
main text. The dashed area is the magnetic ordered phase 
which extends slightly for small Q but split into two parts for 

\n\>2. 



C. Multi-hole Excitations 

It is straightforward to extend the above discussion 
to a multi-hole system. Let us start with p = 2. The 
eigenfunction in the p — 2 sector is a direct product of 
the wavefunctions for the two p — 1 segments, but with 
smaller lattice lengths x — 1 and L — x — 1, respectively, 
here a; = 0:2 — a;i is the distance between the two holes 
sited at xi, X2- 

For the uniform chain, the lowest energy state of the 
p — 2 sector should correspond to the configuration 
where two holes are close together, as in this case the 

total surface energy is minimized. So the lowest hole 

(2) 

excitation gap in the p = 2 sector is given hy A)^' = 
Eq{2,L) — Eo{0,L). For p > 2, one can further show 
that the lowest state in the p-th sector is the configura- 
tion in which all holes are close to each other. This is the 
hole condensation phase in one dimension. The remain 
spin chain has a length L—p. When p << L and L — > 00, 
the finite size effect is negligible so one has A^^^' w pAh- 



D. Phase Diagram 

We now discuss how dimerization influence the phase 
diagram. We have shown that the ground state still lies 
in the p = sector in the presence of dimerization. How- 
ever, dimerization splits the bands and induce a gap at 
higher energies. The quantum critical points extend to 
lines, which are determined by the gapless condition, i.e. 
Ao = 0, or according to Eq. ([?f)) . 



by introducing ( assuming J1J2 > 0) 



R 



D1+D2 
2VJ^ 



D1-D2 
2VJT72 



(30) 



where R describes the competition between Dj and Jj, 
and r2 is a measure of the dimerization strength. In the 
uniform chain limit, R = 2/ A, = 0. With increasing 
il, we have the following three situations: 

(1) In the weak dimerization regime, jfij < 2, there 
is a pair of symmetric critical points, ±i?c, with Rc — 
ViW+A. The ground state is magnetically ordered when 
\R\ < Rc- 

(2) In the strong dimerization regime, > 2, there 
are two pairs of symmetric critical points, ±i?ci,±-Rc2j 
with i?ci = Vri^ -I- 4, i?c2 = Vri^ - 4. The magnetically 
ordered phase appears when Rc^ < R < Rci ■ 

(3) When \n\ = 2 {R^^ = 0), there are three criti- 
cal points, which take values ±2-\/2, respectively. The 
ground state is magnetically disordered when \R\ > 2-\/2, 
but ordered (either ferromagnetic or anti-ferromagnetic, 
depending on the signs of Ji,2) when \R\ < 2V2. 
Note that the R = point corresponds to an al- 
ternating (or staggered) array of single-ion anisotropy 
{D,—D,...,D,—D). It becomes critical when D = 

±277772. 

The R — n ground state phase diagram is plotted in 
Fig. 21 The magnetic ordered phase and paramegnetic 
disordered phase are separated by the critical lines (red). 



VI. THERMODYNAMIC PROPERTIES 



D1D2 = ±4Ji J2 



(29) 



The influence of dimerization can be seen more clearly 



In this section, we study the thermodynamic proper- 
ties of the S = 1 QIM. As was shown previously, this 
model is exactly solvable not only for the ground state. 
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but also for all excited states. However, exactly evalu- 
ating the thermodynamic quantities is still a very hard 
task, particularly for large system size and dimerization. 
Here, we shall develop the recursion method proposed in 
Rcf. 9] in the evaluation of the partition function as well 
as other thermodynamic quantities for either uniform and 
dimerized chains. 



A. Recursion method 

Based on the HDS, the partition function Z{L) of the 
system with lattice length L can be expressed as the sum 
of all the partition functions of the subsystems, i.e., 



ZiL) 



p=0 



(31) 



Where, Z{p, L) is the partition function of the p-th sec- 
tor. Because for fixed p, there are many different hole 
configurations. So Z{p^ L) can be further rewritten as a 
sum over all possible hole configurations 



Z{p,L) 



E Z{xi, 



(32) 



For each hole configuration {xi,--- ,Xp}, the corre- 
sponding partition function of the open chain is given 

by iP = l/ksT) 



Z{xi, ■ 



£'x„)] 



(33) 



Where, H'{ln) is the Hamiltonian of the n-th segment. 

The partition functions of each segments can be 
regarded as the building blocks of total partition 
of the original system. These building blocks are 
denoted by 2:(/„), the partition functions of the 
spin- 1/2 TIM segments with length Then, 
in the uniform case, Z{L) can be expressed as 
Z{L) = Z]p=o ^('i)°^^('2)a---az(/p+i), where 

a = exp(/3D/2) is the partition function of a hole, 
z(l) = 2cosh{(3D/2), and z(0) = 1. As the length of 
the allowed segment may vary, one has the summation 
constraint X]n=i in = L — p. If we denote Z^p'' {L ~ p) = 
z{h)z{l2) • ■ ■ z(Zp+i), then the total partition function can 
be rewritten as 



Z(i)=EafZ(f)(i-p). 

p=0 



(34) 



Where, is contributed from the hole's. 

To numerically evaluate the partition function, it is 
practically convenient to use the following recursion for- 
mula: 



Z(P)(/) = Ez(j-)^(''-i^(? 



-J), 



(35) 



where, Z^-^\l) = 5i,o and Z^°\l) = z{l). By this way, 
we firstly calculate the building block, z{ln), and then 
by iterative use of above relation, evaluate the partition 
function of the S = 1 system. We find that this method 
is particular efhcient for the uniform chain, where the 
system size L could be as large as i = 10000. 

It is non-trivial to extend the above recursion method 
to the spin chain in the presence of dimerization. Here, 
there are four kinds of blocks, associated with the parity 
of the two ends. Thus we can denote them by 2:1.^^2 ('n)j 
with ri_2(= o,e) indicating the left/right ends respec- 
tively. The S* = 1/2 TIM segments with odd or even end 
sites can be solved exactly. The analytical expressions 
for 2;rir2('n) will be provided in a separated supplemen- 
tary material ^285. In the following, we present some nu- 
merical results obtained by the recursion method while 
the system size is kept at L = 2000. For simplicity, we 
use the open boundary condition for the original S — 1 
chain. The extension to the periodic boundary condition 
is straightforward. 



B. The uniform spin chain 

Thermodynamics quantities can be calculated from the 
partition function. In our model, an important physical 
quantity is the thermal average of the hole number, de- 
fined by 



(36) 



Fig. [51^ a) shows the temperature dependence of for 
several A. At low temperatures, Nh increases rapidly 
with increasing temperature for small and positive A. 
The proliferation of the hole number at low temperatures 
is obviously due to the smallness of the hole excitation 
gap. However, for the negative A, say, A = —1, the hole 
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FIG. 5: Temperature dependence of the hole concentration 
Nh (a) and the entropy S (b) for a uniform spin chain with 
A — —1, 0.5, 1, 1.5, respectively. 
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excitation gap is relatively larger, so Nh increases much 
slowly with temperatures. 

Another interesting physical quantity is the entropy 
S, which we plot as a function of temperature for several 
different A in Fig. [Sfb). Wc find that the suppression of 
the entropy S is stronger for larger A. The suppression is 
even more pronounced for the negative A. These behavior 
are similar to the temperature dependence of the hole 
number and are also due to the hole excitation gap. 

We also find that the approaches about 1 /3 in the 
high temperature limit for all A. Correspondingly, the en- 
tropy saturates at the value In 3 in the high temperature 
limit (not fully shown in Fig. [5]). 



C. The dimerized spin chain 
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FIG. 6: (a) Temperature dependence of the hole concentra- 
tion Nh and (b) the entropy S for a dimerized spin chain in 
three representative cases: D1D2 — J1J2 (dashed line), 4Ji J2 
(real line), I6J1J2 (dash-dotted line), respectively. Where, 
Ji = 10, J2 — 1, and Di = 2 are fixed. 

The recursion method, after some extensions discussed 
previously, is also used to evaluate thermodynamic quan- 
tities for a dimerized system. Figs. H] and [7]^ a) show the 
temperature dependence of the entropy S, the average 
hole concentration N^, and the specific heat C in the 
ordered [D1D2 = J1J2), critical {D1D2 = 4J1J2), and 
disordered phases(I?iZ?2 = I6J1J2), respectively. Here 
we fix Ji = 10, J2 — 1, Di = 2, and choose several 
L»2 (= 5, 20, and 80) in Figs. E]and[7);a). Note that 
both the dimerization strength and competition parame- 
ter are tuned by varying D2 ■ The general features of the 
hole number and the entropy are similar to that in the 
uniform case as the hole excitation gap plays the role at 
low temperatures (note that when Ji — 10, J2 = 1, and 
Di = 2, A,, = 7.450, 3.902, and 1.202 for D2 = 5, 20, 
and 80, respectively). 

In Figs. [3 the specific heat is plotted as a function 
of temperature in several cases. We find that at low- 
temperatures the specific heat is peaked at the tem- 
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FIG. 7: (Color Online) (a) The specific heat C for D1D2 = 
J1J2, D2 = 5 (blue-dashed line), L»iL>2 = 4J1J2, D2 = 20 
(red-real line), D1D2 = I6J1J2, D2 = 80 (dash-dotted line) 
with Ji = 10, J2 ~ i, and Di = 2, respectively; (b) The spe- 
cific heat C (dash-dotted line) and the corresponding specific 
heat C(p=0) for the p = subsystem (green open circle) for 
D1D2 = I6J1J2 and Ji = 10, J2 = 1, Di = 2, D2 = 80, re- 
spectively, (c) The specific heat C (blue- real line) and C(p=0) 
(red-doted line) with Ji = J2 ~ 1, Di_ = 2, and D2 — 7. 



perature scale near the hole excitation gap. The peak 
behavior changes depending on the competition of the 
fermionic and hole excitations. The influence of the 
dimerization on the specific heat can be clearly seen at 
relatively higher temperatures, as the dimerization in- 
duced gap is much larger than the hole or fermion ex- 
citation gaps. Fig. Efb) shows the specific heats of the 
S = \ QIM system and its p = sector for the case 
= 10, J2 = 1, Di = 2, £12 = 80. Note that the 
p = sector is identical to the corresponding S" = 1/2 
TIM. In this case Ao('--" Ah) is very small, leading to 
a sharp peak in the 5=1 QIM system. We also find 
a well separated and relatively round peak in the higher 
temperature regime where the energy scale is close to the 
dimerization gap Ad{— 40.025). 

The competitions among the dimerization effect, the 
hole excitations and the fermion excitations can be seen 
more clearly from Fig. ^c) where Ji = J2 = 1, I?i = 2, 
and D2 — 7. In this case, the three kinds of gaps are 
well separated. Ah < Aq < A^, so that the specific heat 
exhibits three peaks. The first one is a sharp peak around 
T ^ Ah [2^, other two peaks are around T ~ Aq and 
T ^ Ad, respectively. By contrast, the low temperature 
sharp peak disappears in the p = sector and other two 
peaks still persist [see the red dotted hue in Fig. Wic)]. 
This is in agreement with the fact that no hole excitation 
exists in the p — sector. 

We also calculated the scaled free energy $s(T) at the 
critical point D1D2 = 4Ji J2 for the dimerized spin chain. 
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FIG. 8: (Color Online) The scaled free energy for the dimer- 
ized systems in the quantum critical phase where D1D2 ~ 
4Ji J2. Two typical situations are plotted: (J2 = 3, D2 — 6) 
and (J2 = 10, D2 — 20), with small and large hole excitation 
gaps, respectively. Other parameters are J2 = 1 and Di = 2. 
As a comparison, the corresponding 5 = 1/2 TIM case is also 
plotted. 



$s(T) is defined by 

^ ^l^-ll^'^'r^'^" (37, 

This quantity, which is identical to the specific heat coef- 
ficient at low temperatures, was introduced in Ref. [sot in 
order to show the temperature persistence of the quan- 
tum critical scaling behavior in the 5=1/2 TIM. In the 
critical region where quantum critical fluctuations dom- 
inate this quantity should be a constant. It was shown 
that in the S = 1/2 uniform TIM it deviates from the 
constant only when T > J/2, indicating a rather higher 
temperature scale below which the quantum critical scal- 
ing behavior persists [s^ . In our recent work Ref . , we 
found that in the uniform 5 = 1 QIM this behavior is 
strongly suppressed by the hole excitations. Here, we 
find that the similar conclusion can be inferred in the 
presence of dimerization. In Fig. [51 we plotted the re- 
sults for two cases (i) Ji = 3, J2 = 1, -Di = 2, £>2 = 6, 
and (ii) Ji = 10, J2 = 1, Di = 2, and D2 = 20. As a 
comparison, the corresponding result for the dimerized 
5 = 1/2 TIM is also plotted. We find that the quantum 
critical scaling behavior at T ^ (or J2/T ^ 00 in the 
Figure [8]) persists approximately at finite T « O.5J2 for 
the dimerized S = 1/2 TIM and the case (ii), and at 
T O.2J2 for the case (i), respectively. This is because 
that the hole gap in the case (i) is much smaller than the 
case (ii). Consequently, the hole excitations play more 
significant role in suppressing the quantum critical scal- 
ing behavior in the former case. This result is consistent 
with the conclusion in Ref. [9| . 



D. Hole condensations: the case with finite Dz 

Now we turn to the case with non-zero D^- It is 
straightforward to show that the ZJ^-term plays the role 
of chemical potential for holes in the 5=1 QIM [§]. 
More specifically, the energy of a p-hole sector with fi- 
nite Dz is related to that with = by the following 
relationship [£^^^'(0) = Eo{p,L)]: 

E^^\d,) = E^^\0)+pD,. (38) 

Thus, with finite D^, the fermion excitation spectra re- 
main unchanged, but the hole excitation gap becomes (3l| 

Ah{D,) = ^h{0) + D,. (39) 
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FIG. 9: (Color Online) Temperature dependence of the hole 
concentration Nu for = or — — Ah(0). Other param- 
eters are Ji = 10, J2 = 1, and Di = 2. 

Therefore, the ground state depends strongly on the 
value of D2. When > — A;i(0), the hole excitation 
is positive, and the ground state is still in the p — 0- 
sector. But when Dz < — A;i(0), the hole excitation gap 
is negative. This indicates that the p = sector is no 
longer the lowest energy state and there are holes in the 
ground state. As a result. Theorem 2 is no longer valid 
in the present case. When L is sufficiently large and p is 
relatively small, the lowest energy of the p-sector can be 

approximated by £^0''^ (0) w i?Q°^ (0)-|-pA;,(0) as discussed 

in the previous section. Then, E^J'^D^) « E^°\Dz) + 
p[A/i(0) -I- Dz], so that we have 

E^o'-HDz) <■■■< E^^\Dz) < E^^\Dz). (40) 

Thus, the order of the band structure hierarchy is com- 
pletely overturned. In this case, the ground state is in the 
p — L sector and all sites are occupied by holes 0. This 
can be also seen clearly from Fig. [9l where the tempera- 
ture dependence of the hole concentration Nh is shown. 
We find that in the zero temperature limit, Nh is equal 
to 1 when < -Ah{0) or when > -Ah{0). 



10 



VII. SUMMARY 

In this paper, we have studied a class of exactly solv- 
able 5 = 1 QIMs with single-ion anisotropy. They exhibit 
a hierarchy of the band structure with both fermionic and 
hole excitations. The hole excitation gap can be tuned by 
the longitudinal crystal field D^- It becomes zero when 
Dz is equal to — A;i(0). The ground state exhibits three 
distinct phases: the magnetically ordered or disordered 
phases when > —Ah{0), or the hole condensation 
phase when Dz < — A/i(0). The transition to the hole 
condensation phase is of the first order. 

We have shown that dimerization does not destroy the 
exact solvability of this model. To our knowledge, this 
is the first example in dimerized 5 = 1 quantum spin 
systems where all the eigen states as well as the wave- 
functions and the eigen energies can be solved exactly. 
The hole excitations enhance the thermodynamic fluctu- 
ations as evidenced in the specific heat which shows a 
sharp peak in the low-temperature region where the hole 
excitations proliferate. This strongly reduces the charac- 
teristic temperature below which the quantum criticality 
persists. All these results are robust against the dimeriza- 
tion. However, dimerization deforms the phase diagram 
and affects the high energy behavior. 



We have developed a recursion method to sum over 
all hole configurations efficiently. This provides a pow- 
erful approach for evaluating rigourously all thermody- 
namic quantities as well as static and dynamic correla- 
tion functions of the QIMs in the thermodynamic limit. 
The detailed derivations for these quantities in the pres- 
ence of dimerization will be provided as a supplementary 
material[28]. The recursion method holds not just for the 
model studied here. It can be easily extended and applied 
in other physical systems whose Hamiltonian can be writ- 
ten as a sum of independent spin segments, separated by 
nonmagnetic impurities, such as Pd- or Zn-doped quasi- 
one-dimensional antiferromagnets Sr2(Cui_a;Pd:r)03 or 
Cui-.Zn.GeOalli,!!!. 
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